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We extend the notion of quasi-exactly solvable (QES) models from po- 
tential ones and differential equations to Bose systems. We obtain conditions 
t^^ . under which algebraization of the part of the spectrum occurs. In some par- 

o, 

On , ticular cases simple exact expressions for several energy levels of an anhar- 

(»— ^ . monic Bose oscillator are obtained explicitly. The corresponding results do 



Dh: 



not exploit perturbation theory and include strong coupling regime. A generic 



^ . Hamiltonian under discussion cannot, in contrast to QES potential models, 

^ . 

^ . be expressed as a polynomial in generators of s/2 algebra. The suggested 

approach is extendable to many-particle Bose systems with interaction. 
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I. INTRODUCTION 

An anharmonic oscillator represents one of the "eternal" problems and models of theo- 
retical physics. It serves as a basis for checking different approximate methods in quantum 
mechanics, the simplified counterpart of field-theoretical models, etc. Apart from this, it 
is of interest on its own since the real world certainly deviates from idealized picture of 
harmonic oscillators due to interaction between them and self-interaction. In so doing, the 
notion "anharmonic oscillator" is mainly applied at least to two different entities. First, it 



refers to some potential power models in which the potential contains terms with the higher 
degrees with respect to coordinate. (The literature on this subject is so vast, that it is even 
rather difficult to indicate some concrete references - let us mention here only the reviews 
|1|], 0, the book and references therein). Second, it is related to quantum Bose models 
with interaction or self-interaction. In both cases the Schrodinger equation cannot be solved 
exactly. However, for the ffist case it was realized that, in spite of impossibility to find the 
whole energy spectrum exactly, in some particular cases (for instance, sextic oscillator with 
special relationship between coefficients [Q) one can find the part of the spectrum (more 
precisely, algebraization of the part of the spectrum occurs). Such a system is the example 
of so-called quasi-exactly solvable (QES) ||5[ ones which includes a rather vast class of poten- 
tials and have direct physical meaning, first of all related to properties of magnetic systems 

i- 

The aim of the present paper is to extend the notion of QES systems to Bose ones and 
apply QES approach to anharmonic Bose oscillators. Strange as it may seem, the approach 
to Bose oscillators in the spirit of QES models was, to the best of our knowledge, absent 
in literature before in spite of the developed apparatus relating realization of Lie algebras 
in Fock space and properties of differential equations |0|. Meanwhile, the QES approach 
to Bose systems deserves treatment on its own due to an obviously wide area of physical 
applications. 

In the coordinate-momentum representation 

a -^ -—, a^ -^ X, (1) 

ax 

Bose Hamiltonian, polynomial with respect to a, a"*" becomes a differential operator and any 
possible invariant subspace is spanned on the polynomial basis. For the systems of such 
a kind there exists the Turbiner's theorem that states that the most general QES can be 
expressed in terms of generators of s/2 algebra 

J+ = a+^a - A^a+, J" = a, J^ = a^^a . (2) 



(the relations (^ are known in the magnetism theory as the Dayson-Maleev representation 
PJ). Nevertheless, the problem of finding QES Bose Hamiltonians cannot be exhausted 
by a simple reference to this theorem since, as we will below, for the typical case under 
discussion, the conditions of validity of this theorem are not fulfilled, so Hamiltonians (except 
some special cases) cannot be expressed in terms of s/2 generators at all. This is the point in 
which Bose QES systems qualitatively differ from potential QES models whose Hamiltonians 
are built with the help of s/2 generators, realized as differential operators |^ . 

Apart from this, even in the cases when Turbiner's theorem does apply to Bose Hamilto- 
nians, it is much more convenient to formulate the conditions of QES - solvability in terms 
of coefficients of an original Bose Hamiltonian directly without resorting to operators J* at 
the intermediate stage. 

II. BASIC FORMULAS 

Consider Hamiltonian 

PO so 

H = Ho + V,Ho = J2 epia^af, V = J2 ^[(a+a)'a^ + (a+)^(a+a)^]. (3) 

p=l s=0 

Throughout the paper we assume that all coefficients of Hamiltonian are real. For Hamilto- 
nian (^ to have a well-defined ground state, one should take po > Sq + 2 independently of 
the relations between coefficients or po = sq + 2 provided Sp^ > 2Asg. In the x-representation 

(S) 

we obtain 

We are interested in the solutions of Schrodinger equation of the type | if)) = Y.n=o ^n I n), 
where | n) is the state with n particles: a~^a \ n) = n \ n). For Hamiltonian (^ subspaces 
with even and odd are not mixed. Therefore, it makes sense to consider them separately. In 
X representation (Q) the wave function of even states $ = J2i=o (^i^h ^i — ^^^- 



It follows from (^ that 

H,<!>i = ai<!>i+i + A$;_i + 7^$;, / = 0,1,..L, (5) 

so 
so 

Pi = J2As2l{2l-l){2l-2y, 



s=0 

po 



7z = E ^P(20"- 
p=i 

We are interested in the possibility of the existence of the invariant basic F2L = 
{1, x'^, x^.-.x^^}. The condition of cut off at / = L reads: ai = 0. 

For odd states the invariant basic F2M+1 = {x, x^,... x^^'^+^}, ^m = a;^™'^^ and 

H^^^ = a^^^+i + (3^^m^i + 7m^m, m = 0, 1, ...M, (6) 

s=so 

s=0 

s=so 

Pm=J2 M'^m + l)2m(2m - l)^ 

s=0 
P=PO 

The subspace with tti < M is invariant with respect to the action of H provided cxm = 0. 

The procedure described above is, in fact, nothing else than the Bose version of quasi- 
exactly solvable (QES) models, applied to an anharmonic oscillator. Now we would like 
to point out why for the case under consideration Turbiner's theorem, in general, does not 
hold, so our formulas cannot considered as particular cases of its realization. The point 
is that Turbiner's theorem implies that the space of all polynomials of a given degree is 
invariant with respect to J*: Fn = {1, x, x'^, ...x^}. Meanwhile, in our case, only subset 
F2N (for even states) or F2M+1 (for odd ones) is invariant, whereas the set F^ is not. Only 
in particular cases, when both conditions a^ = (for even states) and aM = (for odd 
states) are satisfied simultaneously, Hamiltonian does become an algebraic combinations of 



In contrast to |TT|, where differential equations was the object of research, in our paper 



the coordinate- momentum representation (|T]), in which the operator a becomes differential. 



is used as an useful device at an intermediate stage only. In principle, one could rely directly 
on the known formulas of the action of operators a, a"*" on a states with a definite number 
of particles without resorting to the representation (|I]). 

III. EXAMPLES 

Consider Hamiltonain whose off-diagonal part reads 

V = Ao{a^ + a+2) + Ai[{a-^a)a^ + a+^{a+a)] + A2[{a-^a)^a^ + a-^\a+af], sq = 2. (7) 

Now 

ai = Ao + 2lAi + {2iyA2, (8) 

A = 2/(2/ - l)[Ao + (2/ - 2)Ai + (2/ - 2)^2]. 

am = Ao + (2m + l)Ai + (2m + 1)2^2, (9) 

Pm = (2m + l)2m[Ao + (2m - l)Ai + (2m - 1)^^42]. 

First, consider even states. In the simplest nontrivial particular case the invariant subspace 
is two-dimensional, L = 1. Then 

ai = Ao + 2Ai + 4^2 = 0, (10) 

$ = ttQ^o + ai^i and it follows from the Schrodinger equation H^ = E^ that —EaQ+(3iai = 
0, aoao + (71 — E)ai = 0. Taking also into account (|TD|), we obtain: 



E = ^±^^ + 8{A, + 2A2y (11) 

In a similar way, one gets for the three-dimensional subspace (L = 2): 

E^ - (71 + 72)^2 + [7172 - 16(5A? + 52AA2 + 140^2)]^ + 3272(^1 + 4^2)' = 0. (12) 

If Ai = A2 = 0, the low-lying energy levels of an harmonic oscillator are reproduced from 
(|ll|), (p^ . The equation ([l2|) can be solved exactly in the particular case Ai = —4^2: 



E = 0, Ji±yi ± JklL-JiL _|_ 192^2 . For odd states in the simplest nontrivial case M = 1 we 
have 

ai=Ao + 3Ai + 9^2 = 0, (13) 

ao = Aq + Ai + ^2, ^1 = 6ao, 

The conditions a^, = and aAf = are different, so in general the invariant subspace 
exists only for even or only for odd states. However, it may happen that both conditions 
are fulfilled. Thus, for L = 1 = M the compatibility of ([T0| ) and (|13D demands Ai = —5^42, 



Aq = QA2. Then we have the simple explicit solutions for 4 levels of Hamiltonian @: 



^ = 70±7i ± /(TO+Ti)! + 24^2 ^ ± a/^ + 72AI 
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IV. EXPLICIT SOLUTION FOR TWO LEVELS 

If the coefficient a^ = 0, the dimension of the invariant space is L + 1. Meanwhile, it 
may happen that, in addition, /3i_i = 0. Then the two-dimensional subspace spanned on 
$L and $L_i is singled out from the L+ 1 subspace that gives explicit simple exact solutions 
for two levels, whatever large L would be. For Hamiltonian (^ it follows from (||) that in 
this case 

^ = 2L(2L-3)A2, Ai = A2(3-4L), (14) 

aL-i = -2A2, Pl = - 4L(2L - 1)^2, 

E. = ^^^ ± ^ ^'^ -J'-'^' + a,.,p, = ^^^ ± ^I^^^T^ + 8Li2L - 1)^. 
The similar procedure can be repeated for odd states. 

V. GENERALIZATION 

The obvious generalization of QES Bose Hamilatonians arises if Hamiltonian itself does 
not have a "canonical" structure under description but can be reduced to it with the help 

6 



of the transformation HK = KH', where the operator K is some function of a and a"*". 
In particular, it can reahze u — v Bogolubov transformations. In what follows we assume 
that such transformations, if needed, are already performed, so we try to generalize such 
"canonical" forms themselves. Consider the action of Hamiltonian H = Hq + V with Hq 
from d^) and 

^ = E E Asklia+aya"' + {a+)''>{a+ay] (15) 

s=Ofc=l 

on the functions $„ = x^", where g > is an integer, n = 0, 1, ...: 

fco fco 

Ha:^n = 7n^n + E "nfc'^n+fc + E Pnk^n-k, (16) 

fc=l A;=l 

PO 

7n = E^P(^^)^' 
p=l 

so 

otnk = ^AskinqY, 

s=0 
so 

Pnk = E Asknq{nq - l)...{nq - kq + l){nq - kqY, n>k, 
Pnk = 0, n < k. 

Provided the conditions of cut off at ra = A^ are fulfilled, this Hamiltonian can possess 
the invariant subspace F = {$„, n = 0, 1, ...A^}, so the wave function of the system 
$ = J2n=o'^n^n iucludes the states with nq { n = 0, 1, ...N) particles only. The conditions 
under discussion read now as follows. For any given k we must demand a^^i^i^ ^ = 0, z = 1, 
...k, so we have k conditions. As A; = 1,2,.. .A;o, the total number of conditions is equal to 
iT-i = 2 ■ ^^ ^^^ other hand, the number of coefficients Ask is equal to n2 = (sq + l)ko. 
The system can be quasi-exactly solvable if n2 > rii, so 2so > k^. In particular, in accordance 
with examples considered above, one can always adjust the coefficients properly, if /cq = 1. 
The approach considered in the present paper allows extension to many-particle systems. 
In particular, for two pairs of Bose operators (a, a^), (6, b^) one may take 

H = Y.Hl''^h!t\ (17) 

i 

where Hamiltonians Hi and hi are built from operators a,a^ and 6, 6+, correspondingly, and 
have the structure (^ or (plsl) . 



VI. CONCLUDING REMARKS 

During recent years, the class of QES was extended considerably to include two- and 
many-dimensional systems, matrix models , the QES anharmonic oscillator with complex 
potentials ( 0, p. 192; '^^), etc. Meanwhile, it turned out that, apart from these (sometimes 
rather sophisticated and exotic) situations, quasi-exact solvability exists in an everyday life 
around us where anharmonic Bose oscillators can be met at every step. In particular, the 
results obtained can be exploited in solid state or molecular physics, theory of magnetism, 
etc. The approach suggested in the present paper, shows the line along which a lot of second- 
quantized models with algebraization of the part of the spectrum can be constructed. This 
approach can be also extended to systems with interaction of subsystems of different nature 
- in particular, between spin and Bose operator, Bose and Fermi oscillators. 
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